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Using relativistic conformal hydrodynamics coupled to the linear o model we study the 
evolution of matter created in heavy— ion collisions. We focus the study on the influ- 
ence of the dynamics of the chiral fields on the charged- hadron elliptic flow V2 for a 
temperature-independent as well as for a temperature-dependent r//s that is calculated 
from kinetic theory. We find that V2 is not very sensitive to the coupling of chiral fields 
to the hydrodynamic evolution, but the temperature dependence of rj/s plays a much 
bigger role on this observable. 
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1. Introduction 

In any attempt of describing experimental results of ultrarelativistic heavy ion colli- 
sions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider 
(LHC) via a viscous fluid dynamic model, crucial physics input reflecting the prop- 
erties of the flowing nuclear matter must be supplemented, like the equation of state 
(EOS) and transport coefficients such as the shear (j/) and bulk (£) viscosities^. 

The EOS of nuclear matter created at RHIC and LHC can be obtained from 
Lattice QCD simulations and is fairly known in the low chemical potential regime 
of heavy ion collisions. In contrast, the temperature dependence of the transport 
coefficients describing the flow of this matter is not precisely known. Although it is 
expected that for the Quark-Gluon Plasma (QGP) 77 and £ will depend strongly on 
temperature, usually in hydrodynamic simulations temperature-independent values 
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for these coefficients are assumed throughout the entire evolution. The impact of 
the temperature dependence of transport coefficients on momentum ani sotrop ics as 
obtained from fluid dynamic models has been investigated only recently ^ 

The comparison of hydrodynamic models to data would allow, in principle, to 
pinpoint the location of the QCD phase transition or of a crossover from hadronic 
matter to the QGP. In this context, in addition to the hydrodynamic degrees of 
freedom related to energy-momentum conservation, degrees of freedom associated 
with order parameters of broken continuous symmetr ies m ust be considered as well 
since they are all coupled to each other (see e.g. RefsJ^El). 

We s tudy, in th e context of second-order dissipative relativistic fluid dynamics - 
see RefsP^SlSIini f or details, the influence of the long-wavelength modes of chiral 
fields on the expansion of the fireball created in Au+Au collisions at ^snn = 
200 GeV. In particular, we aim at quantifying the effect of the coupled evolution of 
chiral degrees of freedom on the flow asymmetry characterized by t>2, when shear 
finite viscosity is taken into account within a simple microscopic model for the chiral 
condensate. 



2. Conformal fluid dynamics coupled to the linear er model 

In this section we describe the chiral- hydrodynamic model used to c ompu te i>2- The 
quark fluid is described by second-order conformal fluid dynamicJ^^, while the 
chiral dynamics is obtained from the linear a model (LSM). Further details on the 
model and its numerical implementation can be found in RefM 
The classical equations of motion of the LSM are given by 

D,D^ a + ^ = -9Pa (1) 

where <j) = (er, if) and 

Pa = 9 <P a d q f fj * -aMp",*") (2) 

/o(p M , x^) is the Fermi-Dirac function and U is the potential exhibiting chiral sym- 
metry breaking. The quarks act as a thermal bath, leading to an effective potential 

V e (4> a ,T) = Ufa) - dq T J ln[l + e- E ' T ] 



where d q = 24 and E = ^Jp 2 + m? q with m q = g 2 J2 a ( l ) a- 

The classical equations of motion of the chiral fields are solved selfconsistently 
together with the conservation equations for the quark fluid 

De = -(e +p)V A1 ?/ 1 + IPV^ + g(p s Da + p ps ■ Dtt) 

(e+p)Du l = htftdjf- - u l u a d a e) - A^IF^ + ,g(p s VV + p ps ■ W) . 
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where {p s ,p ps ) = p a , D = u^D^ is the comoving time derivative, V M = A^ a D a 
is the spatial gradient, and a^ v is shear tensor. The evolution equation for IF*" 
contains terms up to second order in velocity gradients and reads 

d T n ia = -n ,a V„/ — L™ + — — — (T^ a - (4) 

6U T 1Vli T TirU T 



n <i n «>M if t± Du n d n iQ 



2r 7r 77 2 it T M m t ~ u T J 

where 77 is the shear viscosity and (7V, Ai) are second-order transport coefficients. 
From Eqs. ?? it is seen that the evolution of the chiral fields affects the evolution 
of the quark fluid through the sources in the energy-momentum conservation equa- 
tions; in turn, the quark fluid affects the evolution of the chiral fields through the 
densities p s and p ps , that depend on the local values of the fluid dynamic variables 
T and u". 

In order to compute the temperature dependen ce o f the shear viscosity rj(T) 
in the LSM, we adapt the method described in RefP and employ the linearized 
Boltzmann equation in the relaxation time approximation (see RefP). The shear 
viscosity is given by 

4r f d 3 p p\ 

where r = r(T) is the collision time calculated from the averaged cross sections a 
for quark-quark and quark-antiquark scattering processes including 1 /N c next to 
leading order corrections as 

T 1 = 6/0 (&uu— >uu + &ud^ud + 5\m— >uu + ® U u^rdd 

+ °ud^uci) ■ ( 6 ) 

We refer the reader to 

RefP for details on the calculation of the a's 
The speed of sound (? s and r]/s of the LSM with g = 3.2 is shown in Fig. [TJ 
together with the c 2 s corresponding to Lattice QCLff^ (shown for comparison). The 
value of g = 3.2 is chosen because it leads to a smooth crossover with a drop in 
Cg near the critical temperatur^f] T c which is comparable to the one obtained from 
Lattice QCD simulations. For larger (smaller) values of g, the drop in (? s near T c is 
significantly larger (smaller), as shown in Fig. [21 It is seen that decreasing the value 
of g leads to a softening of the crossover and an increase in T c , while increasing it 
leads to a decrease in T c and a steeper variation of c 2 s with temperature. We note 
that a first order phase transition is obtained if g ~ 3.8, but we will not consider 
this case since Lattice QCD simulationJ^ indicate the QGP-hadron transition is 
a smooth crossover (see Fig. []}. 



a Strictly speaking, there is no T c since we are dealing with a crossover and not with a phase 
transition. However, we will still refer to a T c as an approximate critical temperature around 
which varies rapidly 
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For the second-order transport coefficients we takeP tv = 2(2 - In 2)^ an d 
Ai = 2^ corresponding to the J\f = 4 Super Yang-Mills theory. We use a 13 fm x 
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Fig. 1. Speed of sound of the LSM and that of the model EOS of Laine and Schroder connecting 
a high-order weak-coupling perturbative QCD calculation to a hadron resonance gas through a 
crossover {upper panel); and n/s(T) in the LSM with g = 3.2 (lower panel). 
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Fig. 2. Square of the sound speed of the LSM as a function of temperature, for three values 
of the chiral coupling: g = 3, 3.2, 3.4. 
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13 fm transverse plane, set the impact parameter to b — 7 fm and the initialization 
time to To = 1 fm/c. As initial conditions we use u x = u y = and H xx = H xy = 
jjyy — o. These are typical values for the parameters used as input in hydrodynamic 
simulations of heavy ion collisions at RHIC, and were used in Refs.L^I. The initial 
energy density profile is obtained from Glauber's model, with a temperature Ti = 
330 MeV at the center of the fireball. As a reasonable ansatz we take 7?(ro , r ) = 
and a(ro,r) — /^[l — e~^ r / r °^ ] with tq = 9 fm as initial values for the chiral 
fields. We use the isothermal Cooper-Frye freeze-out prescription with a freeze-out 
temperature Tp = 130 MeV. 

3. Charged hadron elliptic flow 

We will now discuss the results obtained for the charged-hadron elliptic flow V2 
computed from the chiral-hydrodynamic model described above. 

Fig. [3] shows i>2 calculated taking or not taking into account the source terms 
in the hydrodynamic equations, for cither the temperature-dependent computed 
from kinetic theory or a temperature-independent rj/s which is chosen to be the 
averaged value throughout the evolution of the fireball. The results shown were 
obtained with a value g = 3.2 in the LSM, which corresponds to a smooth crossover 
between hadronic matter and the QGP as suggested by Lattice QCD calculations 
(see e.g. RefPl). 

It is seen from Fig.[3]that both for a temperature-dependent and a temperature- 
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Fig. 3. Elliptic flow of charged-hadrons calculated taking or not taking into account the source 
terms in the hydrodynamic equations, for either a temperature-dependent or a temperature- 
independent r)/s (taken as the averaged value throughout the evolution of the fireball). The value 
of the coupling constant is g = 3.2. 



November 28, 2012 2:49 WSPC/INSTRUCTION FILE peralta"krein3 



6 J. Peralta- Ramos and Gastaao Krein 

independent Tj/s, the elliptic flow does not depend strongly on the chiral sources. 
More precisely, we see that for a temperature-independent rj/s, t>2 turns out to be 
practically independent of the dynamics of the chiral sources for px < 3 GeV. In 
contrast, for the temperature-dependent rj/s there are small but visible differences 
in the elliptic flow calculated with or without the chiral sources. Specifically, V2 is 
slightly larger at low px and saturates at a smaller value of px when the chiral fields 
are taken into account as sources for the hydrodynamic variables. 

In order to determine the impact of the value of g on the behavior of i>2 with 
Pt, we also calculated vi for a temperature-dependent rj/s and different values of 
g close to g = 3.2. The results are shown in Fig. 21 It is seen that when the value 
of g increases, v% decreases (this is due to the fact that for larger values of g, we 
find that the average value of rj/s throughout the hydrodynamic evolution becomes 
smaller). Comparing the values of v 2 calculated taking or not taking into account 
the chiral fields as sources for different values of g, it is seen that the influence 
of chiral dynamics on vi becomes smaller as g increases. This supports our earlier 
conclusion that the evolution of the chiral fields have only a small impact on i>2- 

It seems clear that there are difficulties in extracting (an average value of) 
i]/s from data on v 2 - Our results show that uncertainties associated with the de- 
pendence of rj/s on temperature lead to appreciable changes in the curve of «2 
versus px, much larger than the changes stemming from the coupling of fluid dy- 
namic variables to evolving classical chiral fields. The uncertainty associated with 
the temperature-dependence of r//s should be added to the theoretical uncertainty 
that comes e.g. from the initial conditions (for example using Color Glass Conden- 
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Fig. 4. Elliptic flow of charged-hadrons calculated taking or not taking into account the source 
terms in the hydrodynamic equations, for g = 3, 3.2, 3.4 and a temperature-dependent rj/s. 
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sate or Glauber initial conditions) and the freeze-out process, that according to 
recent studies ad d up to an overall uncertainty which can be roughly estimated in 



4. Conclusions 

We have found that the values of vi do not depend strongly on the evolution of 
the chiral fields. Specifically, for a temperature-independent rj/s this dependence is 
negligible for pt < 3 GeV, while for a temperature-dependent rj/s it is appreciable 
but still small even at small pt- ^^^^ 

In line with the results of Refs.^EEl, our results show that not knowing precisely 
the temperature-dependence of rj/s leads to significant uncertainties in attempts 
of extracting this ratio from data on V2 , in addition to the uncertainties that stem 
from the initial conditions and the freeze-out process, among others sources. 

It is worth noting that despite the coupling of chiral sources to the hydrodynamic 
evolution would add further uncertainties to the values of rj/s extracted from data, 
they are not very big and do not affect in a significant way the possible extraction 
of rj/s values from charged hadron elliptic flow data. 

Future work includes improvements such as the consideration of bulk viscosity in 
the hydrodynamic equations and fluctuations of the chiral fiel ds. This latter effect 
would act as noise sources in the classical equations of motioJ ^E3D51 
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